Dark Matter Debris Flows in the Milky Way 
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We show that subhalos falhng into the Milky Way create a flow of tidally-stripped debris particles 
near the Galactic center with characteristic speed behavior, but no spatial features. Using the Via 
Lactea-II N-body simulation, we study the unvirialized component arising from particles that were 
bound in subhalos around the time of reionization but have since been tidally stripped. These 
debris particles constitute a few percent of the local density today and have speeds peaked at 340 
km/s in the solar neighborhood. This spatially- homogeneous velocity substructure has important 
implications for surveys of low-met allicity stars, as well as direct detection experiments sensitive to 
dark matter with large scattering thresholds. 



In the ACDM paradigm, structure formation in the 
early universe begins roughly twenty million years after 
the Big Bang with the gravitational collapse of matter 
around density fluctuations [l]. The first virialized ob- 
jects are Earth-mass halos, which continue to grow in 
size, reaching 10^-10^ Mq around z = 10 [2 . Dark mat- 
ter halos of this size host active star formation at their 
cores. These early stars are the pioneers of future galaxy 
growth for the largest halos, but become sparse relics 
in the less massive halos, which experience a dampening 
in star formation rates as gas is photoevaporated during 
reionization [3-5 . In contrast to their larger counter- 
parts, these small halos remain dark. 

The Milky Way halo forms hierarchically by accreting 
these small, selfbound subunits [6 . While some of these 
halos merge completely and virialize, others do not and 
continue to orbit the Galaxy as satellites. These subhalos 
experience tidal forces leading to substantial mass loss, 
particularly on pericentric passage when the forces are 
most extreme [7l|8]. The stripped particles become a part 
of the Galactic halo and those with shared progenitors 
sometimes retain distinctive phase-space features. The 
stars that formed near the cores of the subhalos may also 
be stripped, leading to spatial or velocity features in the 
stellar halo [l{12j. 

Substructure in star distributions can provide hints 
for potential structure in the dark matter distribution, 
as has already been the case for stellar streams asso- 
ciated with the disruption of dwarf galaxies. Stellar 
streams are spatially aligned and have small velocity dis- 
persions [TO'^ ITT] [T3l4l6 . Their presence strongly sug- 
gests that dark matter streams also formed from material 
stripped from infalling satellites. A dark matter stream 
is dynamically cold, has a one-dimensional morphology, 
and can enhance the modulation fraction of a signal at 
a direct detection experiment pTHlQi . Studies of dark 
matter streams in numerical simulations, however, sug- 
gest that a single dark matter stream that dominates the 
local density is rare [2QH22] . 

In this Letter, we introduce another class of veloc- 
ity substructure that is spatially homogeneous. While 
spatially-homogeneous velocity substructure has been 
discussed with regards to the stellar halo (see [9 for a 



review) , it has not received any attention with regards to 
dark matter. However, it can provide a powerful probe 
for experiments because a generic velocity feature that is 
nearly spherically-distributed will be guaranteed in the 
solar neighborhood. Using the Via Lactea-II (VL2) N- 
body simulation, we will show that a subcomponent of 
the dark matter near Earth is comprised of tidal debris 
with a distinct high-speed behavior, but no characteristic 
spatial features. We refer to this velocity substructure as 
a "debris flow," to distinguish it from a tidal stream. 

Via Lactea-II is a high-resolution N-body simulation 
that follows the evolution of a billion 4.1 x 10^ Mq par- 
ticles from z = 104 to the present [8l|23l[24]. It assumes 
a ACDM universe with the best-fit cosmological parame- 
ters from the WMAP-3 year data release [25 . The simu- 
lation is centered on an isolated halo with no recent ma- 
jor merger that might be a suitable candidate for a Milky 
Way-like galaxy; the host has a mass of 1.94 x 10^^ Mq 
and a tidal radius of 462 kpc at 2: = 0. There are 
20047 subhalos identified at 2; = in the VL2 simula- 
tion with peak circular velocity larger than 4 km/s. Of 
these, nearly all have progenitors at z^e = 9 [26 , which 
we choose as the characteristic time when baryon cool- 
ing in the subhalos becomes inefficient due to the onset 
of reionization. We study tidal debris from this subset of 
subhalos because of its relevance for substructure in the 
old stellar halo; a more complete study of debris from 
the full catalogue of VL2 subhalos is explored in follow- 
up work [27 . 

We identify all particles bound to the subhalo progen- 
itors at Zre and then determine their corresponding po- 
sitions and velocities in the present epoch. A particle is 
labeled as "debris" if it is not bound to a subhalo to- 
day. The top left panel of Fig. [l] is a snapshot at z = of 
the debris from a randomly chosen orbiting satellite. The 
tidally-stripped particles are distributed within a 100 kpc 
radius of the Galactic center. There are relatively larger 
densities of particles near the turning points of the orbits. 
Notice, in particular, that the density is very high within 
~ 10 kpc of the center, where particles are stripped on 
the satellite's pericentric passage. In this region, individ- 
ual streams cannot be distinguished and the overdensity 
is fairly uniform. The top right panel of Fig. [T] shows the 
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particles that have been tidahy-stripped from all the in- 
fahing satehites. The debris shows no spatiahy-locahzed 
structure, but rather fohows the mass-contours of the 
prolate VL2 halo [28 . 

The subhalos at z^e range in mass from 10^ — 10^ Mq, 
and roughly 4000 of these progenitors (approximately 
20%) contribute debris particles today within 45 kpc of 
the Galactic center. Of these, the most massive halos 
contribute the largest fraction. In particular, about half 
of the debris consists of particles from halos with mass 
greater than 5 x 10^ Mq, while 20% comes from halos 
with mass greater than 4 x 10^ Mq at Zre- 

Figure |2] shows the speed distributions for debris parti- 
cles located 30-45 (green), 15-30 (pink), and 5-15 (blue) 
kpc from the Galactic center. For comparison, the solid 
black curve shows the distribution for all VL2 particles 
in a 5-15 kpc shell. Note that these distributions are sep- 
arately normalized to emphasize the differences in their 
shapes. The speed of the debris particles increases for 
radial shells closer to the Galactic center. In the 5-15 
kpc shell, the speeds of the debris particles are peaked at 
^ 340 km/s and the shape of the debris distribution is 
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FIG. 1: (Top left) Distribution of dark matter particles at 
present time that were stripped from an orbiting satellite with 
a mass of 5.6 X 10'^ M© at Zre = 9. (Top right) The distri- 
bution at present time of all tidally-stripped particles. Both 
distributions are projections in the x-z plane, centered on the 
Galactic origin and extending out a radial distance of 100 kpc. 
(Bottom) Fractional density of the debris particles relative to 
the Milky Way halo as a function of radius. The spread in 
density in the 8-9 kpc bin is indicated by the shaded bars (see 
text). 



distinct from that of all the VL2 particles in this shell. 

The radial dependence of the debris speeds is a simple 
consequence of energy conservation. The gravitational 
potential is related to the density distribution of the halo, 
which, assuming an NFW-like profile [29]-[3T], is the fol- 
lowing for the VL2 parent halo: 



p(r) 



(r/r,)7(l + r/r,) 



3-7' 
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where ps = 3.5 x 10 ^ M^pc 



28.1 kpc, and 7 



1.24 are the best-fit parameters for the local density, scale 
radius, and inner radius, respectively |28]. The initial 
velocity of the debris is taken to be the average velocity 
of the parent subhalos at final apocenter. Using these 
assumptions for the initial conditions and the potential, 
we estimate that v ~ 275 km/s at r = 25 kpc and v ~ 370 
km/s at r = 10 kpc; these speeds correspond remarkably 
well with Fig. jSJ 

It is clear, then, that the speed behavior of the debris 
can be explained by the common origin of the particles 
in infalling subhalos. The fact that these particles are 
tidally stripped from satellites also affects the radial, Vr, 
and tangential, Vt^ velocity distributions. The top row of 
Fig. [3] shows the regions of highest density in the Vr — 
Vt plane for the debris particles within 30-45 kpc (left), 
15-30 kpc (middle) and 5-15 kpc (right). There is a 
dramatic increase in the tangential velocity of the debris 
particles from ~ 60 km/s far from the Galactic center 
up to ~ 340 km/s within the inner 15 kpc of the halo. 
At large radii, the radial velocity is centered at = 
and has small dispersion. It is maximized between 20 
and 30 kpc, with a value of ^ 180 km/s. In the inner 
15 kpc, Vr is again centered at 0, but with a much larger 
dispersion. The corresponding distributions for all non- 
debris particles in these radial shells are shown in the 
bottom row for comparison. 

The velocities of the debris particles closest to the 
Galactic center have a large tangential component. This 
is because these particles are tidally-stripped at or near 
the pericenter of their parent halo's orbit at small radii. 
Further from the Galactic center, the velocities of the 
debris flow have a larger radial component because the 
particles are stripped from infalling satellites that have 
not yet reached their turning points. A random selec- 
tion of late-falling high-speed particles in the Milky Way 
should not necessarily exhibit this radial dependence in 
Vr and Vf. Indeed, the behavior observed in Fig. [3] is a 
direct consequence of the tidal origin of the debris. 

To determine whether the dark matter debris has im- 
portant observational consequences, it is necessary to 
characterize both its density and velocity behavior in the 
solar neighborhood. The bottom panel of Fig. [l] shows 
the density of the debris relative to the VL2 halo, as 
a function of radius. From 5-20 kpc, where the speed 
behavior is most distinct from the rest of the halo, the 
debris comprises 1-3% of the Milky Way. This fraction 
increases to 4% at 40 kpc; however, as shown in Fig. [sj 
the speeds of the debris flow are not as high at these 



3 



7 



6 




V (km/s) 



FIG. 2: Speed distribution of the dark matter debris. The 
black hne is the distribution of all VL2 particles in a shell from 
r = 5-15 kpc. The dotted lines are the distributions for debris 
particles located 30-45 (green), 15-30 (pink), and 5-15 (blue) 
kpc from the Galactic center. The dark blue region is the la 
spread about the mean for 100 sample spheres with radius 5 
kpc centered at r = 10 kpc. The light blue region shows the 
maximum and minimum for the same sample spheres. 



distances. From 7.5-9.5 kpc, the debris makes up about 
1.7% of the halo. The fractional contribution of debris in 
this radial shell increases to about 4.4% for particles with 
Earth- frame speeds greater than 500 km/s (in June). 

To determine the local variation in the debris flow den- 
sity, we find the relative density in a hundred sample 
spheres with radius 0.5 kpc centered at 8.5 kpc. The 
dark blue bar in Fig. [l] shows the la spread and the light 
blue extends over the total range of sampled densities. 
Part of this variation is due to the fact that the sam- 
ples are taken in a spherical shell, but the VL2 halo is 
prolate [28]. 

To characterize local deviations in the debris' speed 
behavior, we sample the speed distributions for a hun- 
dred spheres centered at r = 10 kpc with 5 kpc radius. 
The dark blue region in Fig. [2] is the la spread in the dis- 
tributions over these samples, and the light blue region 
shows the minimum and maximum value found in each 
speed bin. The spread in the debris' speed distribution 
has some localized peaks, but its overall shape remains 
remarkably consistent over the entire spherical shell. 

In summary, a subcomponent of the local Milky Way 
halo is characterized by dark matter tidal debris with 
unique speed behavior, but no local spatial structure. We 
introduce the term "debris flow" for this class of spatially- 
homogeneous velocity substructure. Debris flows 

• consist of overlapping sheets, streams, plumes and 
shells created by dark matter tidally stripped from 
inf ailing subhalos. 

• have a spatial distribution indistinguishable from 
the background halo. 



• have a peaked speed distribution and unique radial 
and tangential velocity behavior, as determined by 
the orbital properties of the subhalo progenitors. 

The debris in VL2 that was originally bound in subhalos 
at the time of reionization now constitutes approximately 
a few percent of the local density, and has speeds peaked 
^ 340 km/s in the solar neighborhood. 

The debris flows described in this Letter serve as one 
of the first examples of spatially-homogeneous velocity 
substructure in the Milky Way halo, and should be stud- 
ied in different simulations. For example, the fractional 
density of the debris depends on the number of subhalos 
resolved by VL2 and will vary between simulations with 
different resolutions and initial conditions. Despite these 
variations, the debris flow should be generic to similar 
high-resolution Milky Way simulations like GHalo [32 
and Aquarius [22l |33] , whose subhalo concentrations are 
in good agreement with VL2 [28] . The properties of the 
debris flow will also be affected by the inclusion of the 
Galactic disc. VL2 does not include baryonic physics, 
which will increase the internal velocities of the orbiting 
satellites, as well as the energies of the debris particles. 

It is important to remember the selection bias in the 
subhalos considered in this work. In particular, the parti- 
cles labeled as "debris" originate from subhalos that still 
exist in the present epoch and are bound to the Milky 
Way. The debris does not include particles stripped from 
subhalos that either pass through the Milky Way or are 
completely destroyed during infall. In addition, we do not 
consider particles that were bound at redshifts other than 
Zre = 9. Debris from this redshift provides a good start- 
ing point for understanding debris flows because of its 
relevance for star surveys, as well as direct detection ex- 
periments. However, contributions from other redshifts 
should increase the relative density of dark matter debris 
and is explored in follow-up work [27^ . 

The presence of debris flows can be experimentally 
verified in several different ways. One possibility is to 
look for an imprint of the flow in the local stellar dis- 
tribution. Because the time required for momenta ex- 
change between stars is much longer than the age of the 
Galaxy, the kinematics of old stars encode information 
about their origin. Stars that are tidally-stripped from 
subhalos "trace out" the paths of their dark matter coun- 
terparts [9 . These stars are some of the oldest and most 
metal-poor in the Milky Way because they originate from 
subhalos with inefficient star formation after reionization. 
Simulations of dark matter and baryonic evolution have 
followed the accretion of satellites in simple galactic mod- 
els and have found evidence for a rich morphology of 
structure in the stellar halo [34]. The specific evidence 
for debris flow would be metal-poor stars that exhibit dis- 
tinct velocity behavior, but no distinct spatial features, 
over large areas of the sky. 

There is accumulating experimental evidence for sub- 
structure from surveys that study both the photometric 
and spectroscopic properties of stars over large fields of 
view. One of the most dramatic examples is the discov- 
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FIG. 3: Tangential vs. radial velocity for debris particles within 30 < r < 45 kpc (left), 15 < r < 30 kpc (middle), and 
5 < r < 15 kpc (right) of the Galactic center. The top row shows the distributions for the debris particles, while the bottom 
row shows the distributions for all other particles in VL2. 



ery of the Sagittarius dwarf galaxy [3F and the trail of 
stars that form a stream along its orbit [36-40 . Sev- 
eral other spatial overdensities have been located in the 
stellar halo [41-43 and efforts have been made to detect 
structure in velocity-space j44H46] . For example, a recent 
SEGUE study identified new elements of cold substruc- 
ture from excesses in line-of-sight velocities over stellar 
halo expectations [46]. The upcoming GAIA [47] satel- 
lite will reconstruct the proper motions of stars and will 
be an integral step in mapping the velocity substructure 
in the solar neighborhood. 

Direct detection experiments, which look for nuclei re- 
coiling from a collision with a dark matter particle, pro- 
vide an alternate test for debris flows. In particular, the 
debris results in high-speed flows in the Earth frame that 
yield a distinctive energy spectrum for the recoiling nu- 
cleus [27]. Models in which the dark matter has a mini- 
mum scattering threshold near the velocities of the debris 
particles would be particularly sensitive to this feature in 
the velocity distribution. Light elastic dark matter [48- 
[53j is one example that has recently received considerable 
attention in the literature. 

Debris flows offer a unique way to search for Galactic 



dark matter, with star surveys and direct detection ex- 
periments providing orthogonal detection possibilities. A 
discovery would indicate that a significant fraction of the 
local halo is unvirialized and retains distinctive phase- 
space features, despite being spatially uniform. More 
importantly, it would open a window to the long history 
of our Galaxy's evolution, providing a fossil record of the 
subhalos that merged to form the Milky Way. 
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